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ABSTRACT 

We derive an expression for the correlation function of the random force on 
a sohton which is consistent with the constraints needed to integrate out 
the zero modes which appear due to the broken translational symmetry of 
the sohton solution. It is shown that when the constraint does not com- 
mute with the operator which defines the correlation function, i.e. when 
the operator is not physical, only low frequency phonons contributions 
may be considered. On the contrary, when the correlation function of the 
random force on the soliton is constructed with physical operators one 
may also include in a correct manner the contributions from the optical 
phonons. 
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1 Introduction 



Although the 0^ is not an adequate model to describe solitons and their motion in 
polyacetylene because it does not include retardation effects, it is a model simple 
enough to make possible the study of soliton diffusion In this model, the soliton 
is taken as the classical solution of a non linear differential equation. This solution 
is stable because it is possible to define a topological charge which commutes to 
leading order with the Hamiltonian describing the soliton. It is then impossible 
to deform the soliton solution to a trivial solution in a continuous manner. The 
quantum fluctuations about the soliton are then the phonons of the model. Since the 
Hamiltonian is non linear, there are interactions between the soliton and the phonons. 

One may then consider how these phonons give rise to a Brownian motion of 
the soliton. The observable which measures how the Brownian motion takes place 
is the diffusion operator. The diffusion operator has a diffusion constant which is 
proportional to the square of the phonon number. Thus, at zero temperature we 
expect the soliton to undergo a mild Brownian motion due to quantum effects. As 
temperature increases, soft frequency phonons will contribute to the diffusion of the 
soliton. As temperature increases even more, optical phonons will come into the 
picture and contribute to the Brownian motion of the soliton. Thus, in order to 
describe soliton diffusion at high temperature in a correct manner it is necessary 
to take into account the contribution of optical phonons. Experimental values of 
the extension of the soliton are of the order of the wavelength of optical phonons. 
It is thus important that the information of the profile of the soliton entering the 
correlation function for the random force be in agreement with that of the classical 
solution. 

The problem of soliton diffusion has been studied in [Q, were used was made of 
the formalism presented in |^ to calculate the contribution to the diffusion of the 
soliton due to phonons. However, the analysis done in Q does not take into account 
in a proper manner the shape of the soliton. Thus, the contribution from the optical 
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phonons to the correlation function used in is incorrect. This follows from the 
analysis done below were we show that the correlation function as defined by Mori ||^ 
is constructed from an unphysical operator. This operator is not physical in soliton 
models because these systems have zero modes associated to the loss of translational 
invariance of the classical solution. In order to integrate out this zero mode and obtain 
a theory free of infrared singularities we must introduce a collective coordinate. In 
doing so we enlarge the phase space, and to recover the physical theory we must make 
use of a constraint which appear in the definition of the collective momentum. That 
is, the collective momentum is not the time derivative of the collective coordinate. 
Rather, it is a function of the fluctuations and classical solution. This relation imposes 
a constraint. And this constraint does not commute with the collective coordinate. 
Thus, the collective coordinate is not physical. This can alternatively be understood 
by looking at the original phase space. This space is embedded in a large space which 
also contains the collective coordinate which was not present in the original theory. 
Despite subtleties which involve certain transformations among the variables in the 
enlarged phase space, the collective coordinate is an artificial operator which is not 
physical. 

The reason is that one cannot simply measure the position of the soliton. The best 
one can do is to measure the energy density of the soliton. Thus, one can effectively 
measure the classical soliton solution. Where the solution vanishes we define the 
position of the soliton. In fact, if one expands the soliton solution in a Taylor series 
in the collective coordinate, she or he will find that the leading term is nothing but 
the collective coordinate. However, the soliton solution is not physical either. To be 
physical one must add the fluctuations about the soliton. Thus, one must replace 
the collective coordinates appearing in the correlation function of the random forces 
which act on the soliton defined in by the field appearing in the theory. 

When such a replacement is made, we find that the information entering the 
correlation function of the random forces acting on the soliton will contain all the 



2 



information of the soliton profile and not just the first term in the Taylor series in 
the collective coordinate. Retaining only the first term is equivalent to assuming 
that the energy density of the soliton is evenly distributed throughout the polymer. 
Thus, low frequency phonons will not see the shape of the soliton but neither will 
the optical phonons. While the soft phonon have wavelength much larger than the 
measured extension of the soliton and thus do not effectively see the soliton, this is 
not the case for optical phonons whose wavelength is of the order of the extension of 
the energy density of the soliton. Thus by keeping only the first term in the Taylor 
series, e.i. the collective coordinate, the contribution of the optical phonons to the 
correlation function of the random forces acting on the soliton is incorrect. Only when 
we retain the remaining terms in the expansion will the contribution of these phonons 
be correct. Thus, conclusions regarding soliton diffusion at high temperatures using 
the diffusion equation in 0] must be reconsidered and effects from optical phonons 
recalculated. 



We consider the propagation of a soliton in a box of length L. The Lagrangian is 



and it is invariant under translations. The first, second and third terms on the 
r.h.s refer to the bosonic term, fermionic term and interaction term between them, 
respectively. The fermions will not play a role in our discussion. We thus concentrate 
on the bosonic sector. The Lagrangian for this sector reads 



2 The Model 



L($,^)=L(<l>)+L(^) + L,„t 



(1) 




(2) 



The classical equation of motion for this system is 



(3) 
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where Vn = d'^V/dcf)^. Static solutions to this equations can be classified according 
their translational invariance. For example, the constant solution 

0CO = A (4) 

has translational invariance, and vanishing topological charge. The Solution without 
translational invariance 

X 

(j)c = tanh{-j=), (5) 
satisfies 

0,(x)^0,(x + X)^0(x). (6) 

The solution 0c is stable since it has topological charge equal to one and cannot decay 
without changing this value. 

The loss of translational invariance of the classical solutions lead to the appearance 
of zero modes in the spectrum of the fiuctuations. The spectrum of the fiuctuations 
q with 

<^ = (j), + q. (7) 
is obtained from the linearized equation of motion for the fiuctuations 

-'q + q" + V2U£ = 0. (8) 
Expanding the fluctuations in normal modes 

g = ^— =(V^„a+ + /i.c.), (9) 
we see that among the solutions ipn to there is a normalized solution 



where M is the mass of the soliton. The eigenfrequency of this solution is Ui = 0. 
This is the zero mode. In doing a perturbative treatment of an observable we will 
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encounter infrared divergences coming from the contributions of this zero mode. In 
order to integrate out these zero modes and obtain an infrared divergent free theory 
we must make use of collective coordinates. 
The solutions ipn to the equation (H) satisfy 

s{^-y) = T^M^Wniy) (11) 

n 

Ipn'ipm = ^nm (12) 



With the help of these relations we can write the quadratic Hamiltonian in the form 

H=j:M<ar.+h + ip'^'^r (13) 

n ^ 

where the summation in the first term excludes the zero mode and the last term 
depends only on the zero mode creation and annihilation operators af and ai 

P^'^=^i^{ai + at) (14) 



3 The Collective Coordinates 

We gauge the translational invariance [Q. We do this by raising the parameter X to 
a dynamical variable X{t) which we will refer to as collective coordinate. The field 
will have a dependence on X{t) 

$(x) ^ $(x + X(t)) (15) 

Considering this dependence in the Lagrangian we find that we recover the equation 
of motion for $. In going to the Hamiltonian formalism, the conjugate of $ is 

n := ^ = TT (16) 

where [fc{x), q{y)] = —i6{x — y). We also acquire an additional equation which follows 
from the definition of the collective momentum. 



This is a constraint which allows us to integrate the zero mode without encountering 
infrared divergences. The constraint is best put in the form 

f = j<P'cP + q'^-P (18) 

Standard treatment of systems with collective coordinates requires that physical op- 
erators commute with the constraint and that physical states be annihilated by the 
constraint. The reason for this requirement is simple. We have embedded our theory, 
which was defined in a phase space which contains $ and its conjugate only, into 
a larger phase space which also includes the collective coordinate and its conjugate. 
The embedding is performed by the gauge transformations whose generator is /. We 
require that physical operators, those which were defined in the original theory, be 
independent of the collective coordinate before the embedding. This is equivalent to 
requiring that the operator commute with the constraint after the operator has been 
embedded. 

Examples of operators which are physical are the Hamiltonian H and the field $ 
for which it holds 

[i/,/] = 0, [$,/] = 0. (19) 

Operators which are not physical are the fluctuations q and the collective coordinate 
X since the commutators 

[/,X] = -z [f,q{x)] = -t^'{x) (20) 

do not vanish. To determine if X is a physical operator we first write the Hamiltonian 
after integrating out the zero mode. The infrared free collective Hamiltonian can be 
written as 

Hcoii = -{{P - I q'TT?, ^—f7^}+- (21) 



Using the definition 



we obtain an expression for X 

X = -UP - f q'Tc), Kr^}+. (23) 

It is a simple exercise to check that the commutator 

[f,X]^0. (24) 

Thus, the operator X does not have a physical meaning. A similar exercise for X 
shows that neither this operator is physical. 

4 The Diffusion Equation 

The Fourier transform of the correlation function of the random force on the soliton 
was written in as 

Ti.) ^ [ ^^^m^e--dt. (25) 



Equation (|25|) is related to the diffusion equation for the soliton 

D{^) = . " , V (26) 
Where a is the thermal average of the square of the soliton velocity. For adiabatic 



phonons, equation (|25| ) is a good approximation. The reason for this follows from the 
fact that the operator X and its time derivatives are not physical operator. Thus the 
equations (|25|) and (^6]) do not have a physical meaning. 



The natural question which arises is which is the minimal operator that must be 
added to the operator X to have a physical operator. To leading order, the constraint 
implies 



P 



J (27) 



Since P is the conjugate of X it follows that 

m = (28) 



The operator in the r.h.s. of p8| ) is unphysical and the minimal operator which 
includes the r.h.s of (123) and is physical at the same time is 



(j) = (pc + Qre + i'lX{t) = 0c + Qre + J = (pc + Q- (29) 

where in the second term use was made of (0) and 

Qreix) = ^iJn{x) / (30) 

The equation ( ^9]) is a sensible choice because the equation (|26|) refers to the diffu- 
sion of the soliton and not of the collective coordinate of the soliton. However, the 
collective coordinate is a good but unphysical approximation to the soliton. To see 
that this is an approximation we write the equation (^) in terms of the soliton field 

n.)=I^Mm^e-^dt. (31) 
^ ' J < $(t)<l>(0) > ^ ^ 

When we associate X with $ we have that 

X = $ (32) 
which to leading order yields after use of (^) 

X = (33) 
The equality holds exactly when 

(pcix) = X (34) 

which is the expansion of the soliton field about the kink position. Thus when we use 
( pSD we are taking (0) to be the solution of the field at arbitrary distances of the 
kink position, rather than the solution (^. 
Using the identity 

< P'\(j)c{x + X)\P>= J dz(j),{z - x)e'^^'^e'^^^'* = (f),{AP')e'^''p'' (35) 
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where AP' = P — P' and AEp/ = ^ ? ^^^^ equation ( pTj) becomes for 

the forward term 

r r.;^ - (df IdP'AE%\MAPTe'^^-' .... 



. ' J dP"AEl„\(j)^{AP")\'^e^^^p"' 
The backwards component reads 



It follows from the construction of T{w) that the equation ( Pq) for the soliton diffusion 
is mediated by phonons whose wavelength is AP 0. 

We know that soft acoustic phonons will not see the soliton profile because their 
wavelength is much larger then the soliton extension. On the other hand, we see 
that the equations (|37|) and ([36| ) depend on the soliton profile. Thus if we replace 
(pc by X we must impose a cutoff on the integrals over the momenta P' P". The 
cutoff should be such that the smallest wavelength is much larger than the soliton 
extension. However, if we do not make such a replacement, we may safely integrate 
over the phonon spectrum, optical phonons included. 

To see the difference between using the classical solution and using the collec- 
tive coordinate in the diffusion equation (pB]) it is convenient to write ( ]26|) as 0] 



Dp{uj) = J dt < P\X{t)X{0)\P > e-"^\ (38) 

If we replace the collective coordinate of the soliton with the soliton field in ( ^8]) we 
obtain 

Dp{uj) = J dt < P\^t)^{0)\P > e'"^'. (39) 

which after using the integral representation for the Heavy-side function and inte- 
grating in time takes the form 



1 ^^AEhm^n 
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Dp{u) = - dP' P'^l^' . (40) 

27r J^Pc AEp, — uj^ 



We may then study how equation ( ^Of ) behaves for different profiles as a function of 
the cutoff. First we note that since we are interested in the contribution of optical 



phonons and as pointed out in |[^, only the leading order term in the uj expansion 
represent the contribution from the optical phonons, we may discard terms which 
depend on uo. Unfortunately, we cannot get an analytical expression for the profile 
(^. However, assuming that P = 0, we may study how -Do(O) = D behaves as a 
function of the cutoff, Pc, in the integral (^). The behaviour is plotted in Figure |1|. 



D 




Figure 1: Arbitrary Units 

We find that the diffusion equation saturates for high values of Pc and it is inde- 
pendent of Pc for large Pc allowing us to carry the integration of D to arbitrary large 

Next we investigate the behaviour for the profile $ = X. It will give us a taste of 
how sensible is the diffusion equation (|39| ) to the profile of the soliton. In this case 
the contribution from the optical phonons is obtained analyticaly 

D ~ dP'sin^{-LP') (41) 
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where L >> 1 is the length of the box in which the sohton is confined. Evaluation of 
the integral as a function of the cutoff shows that D grows like 

Dr^P, (42) 

This dependence was expected and shows that an improper treatment of the collective 
coordinates leads to a diffusion equation which is ill defined. We also conclude that 
D is quite sensible to different profiles which lead to different contribution from the 
optical phonons. For the soliton profile we see that the contribution of optical phonons 
saturates for large momenta. However for the profile X, D continues to increase 
linearly in P^- Define as cutoff the value of momenta for for which D saturates when 
considering the correct soliton profile. Then this must be the cutoff imposed on the 
profile X to have a sensible theory. Thus for acoustic phonons the diffusion equation 
is not sensible to the profile. However, for optical phonons the diffusion equation is 
very sensitive to the profile. 

5 Conclusions 

The existence of a constraint in the system we have considered is a consequence of the 
broken translational symmetry of the vacuum, e.i. by our choice of classical solution. 
It has been shown that the operators defining the diffusion equation for the soliton 
must commute with the constraint. Because the operators defining the diffusion 
equation in are not physical because they do not commute with the constraint it 
follows that this diffusion equation does not have a physical meaning. In order to 
have only the contributions from the physical sector of the theory we must replace the 
operators used in 0] to define the soliton diffusion by physical ones. The canonical 
operator which is physical and contains the information of the soliton position is the 
soliton field. We have also shown that using the collective coordinate as the operator 
defining the soliton profile is not only unphysical but it is also a good approximation 
to the soliton profile if we consider the contribution of the acoustic phonons only. On 
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the other hand, the contribution from the optical phonons, as we have shown, is very 
sensitive to the sohton profile. One must therefore use the exact profile of the soliton 
and not just the linear term in the expansion of the soliton profile when considering 
the contribution of the optical phonons to the diffusion of the soliton. 
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